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ABSTRACT 

We present high quality spectroscopic data for two massive stars in the OB 10 associ- 
ation of M31, OB 10-64 (B0 la) and OB10-WR1 (WC6). Medium resolution spectra 
of both stars were obtained using the ISIS spectrograph on the William Hershel Tele- 
scope. This is supplemented with HST-STIS UV spectroscopy and Keck I HIRES data 
for OB 10-64. A non-LTE model atmosphere and abundance analysis for OB 10-64 is 
presented indicating that this star has similar photospheric CNO, Mg and Si abun- 
dances as solar neighbourhood massive stars. A wind analysis of this early B-type 
supergiant reveals a mass-loss rate of M = 1.6 x lO _6 M0yr _1 , and v^ = 1650 kms -1 . 
The corresponding wind momentum is in good agreement with the wind momentum 
- luminosity relationship found for Galactic early B supergiants. 

Observations of OB 10-WR1 are analysed using a non-LTE, line-blanketed code, 
to reveal approximate stellar parameters of \ogL/L Q ~5.7, T* ^75 kK, v^ ~ 3000 
kms -1 , M/(M Q yr _1 ) ~ 10~ 43 adopting a clumped wind with a filling factor of 10%. 
Quantitative comparisons are made with the Galactic WC6 star HD 92809 (WR23) 
revealing that OB 10-WR1 is 0.4 dex more luminous, though it has a much lower 
C/Hc ratio (~0.1 versus 0.3 for HD 92809). Our study represents the first detailed, 
chemical model atmosphere analysis for either a B-type supergiant or a WR star in 
Andromeda, and shows the potential of how such studies can provide new information 
on the chemical evolution of galaxies and the evolution of massive stars in the local 
Universe. 
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1 INTRODUCTION 

The observation and analysis of hot, luminous stars in the 
Milky Way and other Local Group galaxies allows the in- 
vestigation of stellar evolution and mass-loss within differ- 
ent environments. For example, Wolf-Rayet (WR) stars, the 
chemically evolved descendents of massive OB stars, provide 



* On leave of absence at the Isaac Newton Group of Telescopes, 
Apartado de Correos 368, E-38700, Santa Cruz de La Palma 



keys to our understanding of massive stellar evolu tion and 
nucleosynthetic processes, e.g. Dessart et al. (2000). 



Studies of the physical nature of the A and B-type su- 
pergiants can also place constraints on evolut ion m odels in 
the upper main- sequenc e regi ons, e.g. Venn (1995), McEr- 
lean, Lennon & Dufton (1999). A further use of the A and 



B-type supergiants is in using photospheric analysis to con- 
strain the chemical composition of the interstellar medium 
in their host galaxies. They are visualy the brightest quasi- 
stable objects in galaxies, and being descendants of mas- 
sive OB-type main-sequence objects have small lifetimes (5- 
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20 Myr) . Many elements visible in their photospheres should 
not be affected by internal mixing and contamination with 
core processed material (e.g. O (only slightly affected in AB- 
types), Mg, Si, S, Al, Ca, Fe, Ti, Cr, Ni; Venn et al. 2000, 
Smartt, Dufton & Lennon 1997). Further they can also pro- 
vide new information on the distances of galaxies through 
the detailed study of the strengths of their radiatively driven 
winds using the wind momentum - luminosity relationship 
(WLR, e.g. Kudritzki et al. 1999). An investigation of the 
stellar wind properties of 0-, B- and A-supergiants in Local 
Group galaxies with well defined distances will, therefore, al- 
low us to test the concept of the WLR and its applicability 
for distance determinations. 

Numerous surveys have identified OB and WR st ars be- 
yond the Ma gellan ic Clouds, e.g. Massey et al. (1986); Mof- 
fat & Shara (1987), although little quantitative analysis has 



discussed separately below in Sect, 
general discussion presented in Sect. 
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been carried out to date. The only detailed studies of Wolf- 
Rayet stars beyond the Magellanic Clou ds ha ve been studies 
of late WN stars in M33 by S mith et al. fll997| ) and Crowther 
et al. ( 1997 ). Bianchi et al. ( |l994J ) have published UV spec- 
troscopy of OB-typ e supergian ts in M31 and M33, while 
Monteverde et al. (1997; 2000) derived the oxygen abun- 
dance gradient in M33 from studies of B-type supergiants. 
An indication of the capabilities of 8 m cla ss telescopes has 
been presented by Muschielok et al. (1999) for three B-type 
s uperg iants in NGC 6822 using VLT. Recently, Venn et al. 
(2000) have analysed four A-type supergiants in M31, the 
first detailed stellar analysis in this galaxy. 

The present work demonstrates the capabilities of 4m 
class telescopes by analysing the spectra of two massive stars 
in the OB 10 association of our nearest giant spiral, An- 
dromeda (= M31 = NGC 224). OB 10 (van den Bergh 1964, 
Massey et al. 1986) is located at 23.6' from the centre of 
M31, and has an apparent size of 2x1.0'. Assuming a dis- 
tance to M31 of 783 kpc (Holland 1998, Stanek & Garnavich 
1998, assumed throughout this paper) the association's pro- 
jected size is hence ~ 450 x 225 pc . Its stellar content was 
investigated by Massey et al. (1986) (and later Massey et al. 
1995), who found two Wolf-Rayet candidates, OB 10-WR1 
(classified as WC6-7) and OB 10-WR2 (classified only as 
WN) and a late O-type supergiant (OB 10-150; 08.5Ia(f)). 
The presence of such massive stars implies that it is a very 
young association, having formed in the last ~5 Myrs. The 
de-projected position of OB 10 would suggest it has a galac- 
t ocent ric distance of 5.9 kpc. Blair, Kirshner & Chevalier 
( 1982 ) have determined an abundance gradient in M31 from 
a combination of observations of H n regions and supernova 
remnants, deriving —0.03 ± 0.01 dexkpc" 1 . Hence, OB 10 is 
expected to have a metallicity of ~9.1dex, a factor of two 
greater than the Galactic solar neighbourhood value. 

In this paper, William Hershel Telescope (WHT) spec- 
troscopy of OB 10-64 (B0 la) and OB 10-WR1 (WC6) are 
combined to reveal the actual metal abundance of this re- 
gion of M31, plus the properties of WC-type stars in a galaxy 
simi lar to our own. Designations are taken from Massey et 
al. (1986); OB 10-64 has an a lterna tive designation of 41- 
2265 from Berkhiujsen et al. ( [l988| ), whil e OB 10-WR1 is 
also known as IT5-19 by Moffat & Shara ( |l987J ). 

Our observations are presented in Sect. El Since the 
methods used to analyse the spectra of the M31 B-type su- 
pergiant and Wolf-Rayet star differ considerably, they are 



2 OBSERVATIONAL DATA AND REDUCTION 

WHT spectroscopy was obtained for OB 10-64 and OB 10- 
WR1 during Sep-Nov 1998 as part of a programme to ob- 
serve early-type stars in the spiral galaxies M31 and M33. 
For OB 10-64 these data were supplemented by Keck I high 
resolution echelle spectrograph (HIRES) data around the 
Ha region of the spectrum and Hubble Space Telescope 
(HST) Space Telescope Imaging Spectrograph (STIS) UV 
spectra. 

HD 167264, a Galactic early B-type supergiant, was fur- 
ther observed in September 1998 with exactly the same 
instrumental setup. For comparison with OB 10-WR1, 
HD 92809 (WC6) was observed at the Mt Stromlo & Sid- 
ing Spring Observatory (MSSSO) in Dec 1997. Table|l| sum- 
marizes the M31 spectroscopic observations and published 
photometry. 



2.1 WHT spectroscopy 

The double armed spectrometer ISIS spectroscopy was used 
to observe OB 10-64 and OB 10-WR1 on the nights of the 
28-29 September 1998. The spectrograph slit was positioned 
so that it included both stars, their separations on the sky 
being only 25". Only the blue arm of ISIS was used for the 
first night, with the R1200B grating and a thinned EEV42- 
80 CCD (with format 4096x2048 13.5/xm pixels). On the 
second night, a beam splitting dichroic was used to feed both 
the blue and red channels; for the latter, the grating was the 
R1200R with a TEK CCD used as a detector. For the blue 
spectra, the unvignetted wavelength range was from approx- 
imately 3985A to 4725A with a pixel size of 0.22A; the cor- 
responding values for the red spectra were 5650-6400A and 
0.4A. 

Since the above setup was primarily geared towards use- 
ful photospheric lines for the B-type supergiant, we missed 
the crucial WC abundance lines between 5000-5600A. Sub- 
sequently, service observations were taken of OB 10-WR1 on 
November 24 th 1998, at a lower but sufficient spectral reso- 
lution. The R300B and R300R gratings were used in the blue 
and red arms at dispersions of 0.9A pix -1 and 1.6A pix -1 
covering 3640-6227A and 6404-9277A, respectively. 

The CCD frames were reduced to wavelength calibrated 
spectra using the IRAF reduction systemFI Standard pro- 
cedures were used to bias correct and flat field the stellar 
images. There is considerable background Balmer line emis- 
sion from the nebular region of OB 10, which can vary over 
small spatial scales along the slit. Although nebular emis- 
sion at the sub-pixel level contaminates our stellar Balmer 
line observations of OB 10-64, the CCD images showed that 
the H<5 and H7 regions were less affected than the Ha line. 



T IRAF is written and supported by the IRAF programming 
group at th e National Optical A stronomy Observatories (NOAO) 
in Tucson ( http://iraf.noao. edu|). 
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Table 1. Observing log and fundamental properties of the M31-OB 10 targets. UBV photometry is taken from Massey et al. (1986), and 
My for each star is calculated assuming a distance modulus 24.47±0.07 (Holland 1998, Stanek & Garnavich 1998). The spectral range 
quoted for the ISIS data is for the unvignetted region; the full blue range was not usable due to the large size of the EEV42-80 CCD 
and the small entrance window to the ISIS cameras. The resolution is estimated from the width of arc lines measured at the adopted 
slit-width. The signal-noise-ratio (S/N) per pixel is measured in the continuum after binning by 2 pixels in the spectral direction for the 
ISIS blue arm, and unbinned data in the red. Note that only 1.0 hour integration was obtained for the red ISIS arm in Sept 1998. The 
S/N per pixel in the Keck HIRES data is after binning to 0.3Apix — 1 . 



Star 



Epoch Spectral Ranges (Res) S/N Exp. V 

A Hours mag 



B-V U-B E(B-V) M v 

mag mag mag mag 



WHT ISIS data: 

OB 10-64 Sep 1998 

OB 10-WR1 Sep 1998 

OB 10-WR1 Nov 1998 

Keck I HIRES data: 

OB 10-64 Oct 1999 

HST STIS data: 

OB 10-64 Oct 1999 
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5.5 
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6404-9277 (7) 5-15 
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6530-6630 (0.2) 


~60 


1 






1120-1716 (1.5) 


~20 


2.3 







-0.97 0.18±0.08 -6.93±0.4 
-0.45 0.18±0.08 -5.4±0.4 



Table 2. Equivalent widths, W\, for selected lines of He II, Si III 
and Si IV in OB 10-64 and HD 167264. These lines are very sensi- 
tive to T e g , and the similarity of the line strengths indicate that 
the effective temp erat ures of the two stars are similar to within 
500K (see Section 3.1). Note that the Siiv line at 4089 is blended 
with a nearby O II feature, but at this T B g the Si IV contribution 
dominates and the ratio of the line strengths is approximately 4:1. 
The equivalent widths for OB 10-64 should normally be accurate 
to 10%; and for HD167264 accuracies are normally 5%. 
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225 


Sim 


4567 
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Sim 


4574 
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Si iv 


4089 
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Si iv 


4116 




100 


135 


Si iv 


4654.14 




SO 


80 



Further, our spectral resolution is high enough that the im- 
portant wings of H<5 and H7 will be unaffected by the back- 
ground nebulosity. However, there may be residual errors 
in the cores of the stellar Balmer line profiles. The neb- 
ular emission around WR1 is particularly intense and ac- 
curate subtraction proved very difficult, however as shown 
in FigureH most of the broader WR features can be mod- 
elled while recognising the narrower nebular contamination. 
The spectra were wavelength calibrated using CuNe+CuAr 
lamp exposures that interleaved the stellar spectra. Individ- 
ual stellar exposures were then cross-correlated to search 
for significant wavelength shifts - none were identified. The 
spectra were then combined using SCOMBINE and a variety of 
rejection criteria; the different methods used to combine the 
data were found to have little effect on the signal-to-noise 
and the CR rejection success-rate in the final spectra. 

Unfortunately conditions in November 1998 prevented 
an absolute flux calibration for OB 10-WR1. Nevertheless, 
a relative calibration was achieved, using the spectrophoto- 
metric DO white dwarf standard star Feige 110. This star 



was observed for 300 sec with an identical setup, immedi- 
ately prior to OB 10-WR1 and at a similar air mass. 



2.2 Keck I HIRES and HST STIS spectroscopy of 
OB 10-64 

OB 10-64 was observed with the 10m Keck I telescope, using 
HIRES on 27th September 1997 and again on the 6/7 Octo- 
ber 1999 giving a total of 4 x 3600s exposures. A 1.1" slit was 
used giving a resolution of approximately 35 000, yielding a 
S/N of ~60 at Hq after co-addition and binning in the spec- 
tral direction to 0.3Apix _1 . The CCD echelle spectra were 
reduced in the standard mann er, us ing the same methods as 
described in McCarthy et al. (1997). The most difficult task 



was removing the strong background Ha nebular emission 
from the stellar spectrum. In fact complete removal proved 
impossible, and we have left in the residual (negative) neb- 
ular Ha line in Figs|6| and []. This led to an unrecoverable 
0.9A region of the spectrum, which is however not critical 
to our fitting of the much broader stellar Ha profile. 

HST spectra were taken with the STIS, as part of an 
extensive Programme (G07481) to observe the UV wind- 
lines of massive blue supergiants in M31 and M33. Two ex- 
posures giving a total of 8500s were taken with the FUV- 
MAMA detector on 1st October 1999. The G140L grating 
and the 52x0.2 arcsec slit aperture were used. The data 
were pipelined processed, (including wavelength and flux 
calibration) by the STScI On-The-Fly calibration systemQ 
This processes the data with the most up-to-date calibra- 
tion files, parameters, and software. Further details of the 
observational material can be found in Tablehl. 



2.3 MSSSO spectroscopy 

The Galactic WC6 star HD 92809 was observed with the 
Double Beam Spectrograph (DBS) at the 2.3m MSSSO on 
1997 Dec 24-27. Use of a suitable dichroic permitted si- 
multaneous blue and red observations of HD 92809 cover- 
ing 3620-6085A (300B) and 6410-8770A (316R), plus 3240- 



•f see tittp: //www. stsci.edu/instruments/stif 
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Figure 1. Selected spectral regions for OB 10-64 (solid line) and HD 167264 (dotted line), showing temperature sensitive features (Hell 
lines at 4541 and 4686A; Si III lines at 4553, 4568 and 4574A and Si IV lines at 4089 and 4116A) and the surface gravity diagnostics 
(Balmer lines: H54101 and H74340A). The high quality data of this 18 m star allow us to identify and measure relatively weak stellar 
absorption lines due to metals and He I. 



4480A (600B) and 8640-1 1010A (316R). A 2" slit and 
1752x532 pixel SITe CCDs provided a 2 pixel spectral res- 
olution of ~5A. Wide slit observations were also obtained 
for HD 92809 and the HST spectrophotometric standard star 
HD 60753 (B3 IV) to achieve reliable flux calibration. Atmo- 
spheric correction was achieved by observing HR4074 within 
an air mass of 0.10 from HD 92809. A standard data reduc- 
tion was again carried out with iraf. 

These observations were supplemented by high reso- 
lution (HIRES), International Ultraviolet Explorer (IUE), 
large aperture, short-(SWP) and long-wavelength (LWP, 
LWR) ultraviolet observations of HD 92809. These were 
combined to provide a single high S/N dataset, flux cal- 
ibrated using the curve obtained by Howarth & Philips 
(1986). Finally, the spectra were tra nsferred to the starl ink 
spectrum analysis program DIPSO (Howarth et al. 1996 ) for 
subsequent analysis. 



2.4 Photometry, distance and reddening to M31 

The published photometry of both stars (from Massey et al. 
1986) is provided in Tab le [j]. Assuming (B-V)„ = -0.26 
for a B0 la atmosphere (Deutschman et al. 1976), we es- 
timate an interstellar reddening of E(B— V)=0.18 mag to- 



wards OB 10-64. Massey et al. (1995) have estimated the 
mean reddening toward the massive stars with spectroscop- 
ically confirmed types in OB 10 and OB8 + OB9, deriving 
a value of E(B-V)=0.16 ±0.02. The reddening directly to- 
wards OB 10-64 is hence compatible with the me an va lue 
for the region. Assuming the distance of Holland (1998) to 
M31 of 783±30kpc implies M v = -6.93 ±0.42 for OB 10-64 
(from the errors quoted in TablellJ). We discuss the impli- 
cations of this in terms of wind and evolutionary models in 
Section H. 

Our HST STIS data was combined with existing archive 
FOS and GHRS data to give a flux distribution from 1120- 
3300A. For the value of E(B-V)=0.18, the de-reddened 
HST/FOS flux distribution for OB 10-64 matches theoret- 
ical models for A > 2000 A but shows a noticeable drop at 
shorter wavelengths, possibly indica ting a differen t far-UV 
reddening law to the Galactic case (Seaton 1979). Bianchi 



et al. (1996a) investigated the extinction in M31 using the 
HST/FOS data for OB 10-64 and additional M31 OB stars 
and concluded that the 2175A bump is weak or absent in 
M31. The slope of the extinction curve was consistent with 
the Galactic curve, but the low reddening of their M31 tar- 
get stars prevented strong constraints being placed on the 
exact shape of this extinction curve. We have briefly re- 
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investigated this issue using a slightly different procedure, 
following Fitzpatrick (1986) in using other B-type super- 
giants in the LMC as the comparison stars. We chose Sk- 
67108 (BO la) and Sk-69 276 (09-BOIa) as the comparison 
stars on the basis of the morphology of their UV spectra 
compared to OB 10-64. These two stars have essentially lit- 
tle or no internal LMC extinction and the E(B-V) values 
(0.06 and 0.10 respectively) reflect the reddening from in- 
tervening Galactic material. This Milky Way extinction is 
similar to what we expect for OB 10-64. Using these two 
stars and OB 10-64 we derive a mean M31 extinction law, 
which confirms the lack of a 2175A feature in M31, similar 
to the SMC (Prevot et al. 1984). We find a slope for M31 
which is marginally shallower than the SMC, but steeper 
than that of the Milky Way. A linear fit of E(A-V)/E(B- 
V) versus A -1 (/i _1 m) gives a slope of approximately 1.8 
for M31; compared with approximately 2.2 for the SMC. 
This analysis suggests that the far-UV extinction in M31 is 
greater than that predicted by the Galactic extinction law, 
and that there is a lack of the 2175A bump. However this 
should be treated with some caution as we have only one 
star available here, and it is only lightly reddened by dust 
in M31 itself. Ideally stars with higher M31 extinction (or 
preferably a range in values) are needed to investigate the 
relationship fully. 

For the Wolf-Rayet star, broad-band Johnson photom- 
etry are (differentially) contaminated by strong emission 
lines. Consequently, the Massey et al. (1986) V-band magni- 
tude for OB 10-WR1 will overestimate the true continuum 
flux, and the observed B-V colour does not permit the de- 
termination of interstellar reddening. Therefore, we assume 
an identical reddening to that derived above for OB 10-64 
and estimate a Smith narrow u-band magnitude for OB 10- 
WR1 as follows. We have convolved our WHT spectra of 
OB 10-WR1 and OB 10-64 with a Johnson V-band filter, 
revealing AV=1.25 mag, in excellent agreement with 1.22 
from Massey et al. (19861). Applying the «-band filter to the 
spectrum of OB 10-WR1 yields (v— V)=0.35 mag, so we es- 
timate u=19.67 mag. This agrees well with the narro w-band 
4752A continuum magnitude of 19.6 quoted by (Massey 
& Johr json 199^ ) for OB 10-WR1. Taking the above dis- 
tance and adopting a reddening of E(B — V) — 0.18 implies 
M„ = — 5.4±0.4 mag. This result is significantly different 
from the mean absolute magnitudes of M v = — 3.7 mag and 
—4.8 mag for respectively, Galactic WC6 and WC7 stars 
(e.g. Smith et al. 1990). 



2.5 Spectral types 

While spectral typing is often neglected when carrying out 
quantitative atmospheric and wind analysis of massive stars, 
we quote our revised types to allow the stars to be put into 
contex t with other morp holog ical work. Humphreys et al. 
(1990) and Massey et al. (1993) gave a spectral type of Bl I 
for OB 10-64 based on spectral data of 2— 3.5 A resol ution 
and moderate signal-noise, while Bianchi et al. (1996b) pre- 



confirms the WC6 classificat ion u sing either the Smith et 
al. ( 1990j ) or Crowther et al. ( 1998) schemes, plus the strong 
nebular contamination. 



3 B-TYPE SUPERGIANT ANALYSIS 

The non-LTE model atmosphere techniques have been de- 
scribed in detail by McErlean et al. (1999 - hereafter des- 
ignated as MLD) and here we only provide a summary. 
A grid of non-LTE m odel atmosphe res was generated us- 
ing the code tlusty (Hubeny 1988) for effective temper- 
atures, T c ff, upto 35 000 K and logarithmic gravities from 
log g = 4.5 down to near the Eddington stability limit. Mod- 
els were calculated for two helium fractions, viz. y — 0.09 
(solar) and y = 0.20, where y=N[Re]/N[¥L+Re\. The models 
omit a number of physical processes, including metal line- 
blanketing and wind effects. Furthermore, the assumption 
of a plane-parallel geometry may be of limited validity for 
the low gravity objects considered here. The consequences 
of these omissions have been discussed in MLD. 

The line forma tion calculation s were performe d using 
the c odes DETAIL ( |Giddings 198l| ) and SURFACE ( |Butlei 
1984). Microturbulent velocities, which are close to the speed 
of sound, have previously been found for B-type supergiants. 
Therefore, in the calculation of line profiles, microturbulence 
has been included as an extra free parameter. Metal ion pop- 
ulations and line-profiles were calculated using mainly the 
atomic data of Becker & Butler (1988, 1989, 1990). Such cal- 
culations explicitly include the effects of the relevant ions on 
the radiation field. However they do not include the effects 
of line and continuum blanketing from other metals. 

Significant difficulties were encountered in running DE- 
TAIL particularly for the silicon model ion but also for 
some other species. These difficulties occurred mainly at the 
lowest gravities for effective temperatures greater than ~ 
25 000 K. Examination of the line profiles indicated that this 
was caused by either emission or 'filling in' of the profiles. 
Indeed the DETAIL calculations normally showed an over- 
population of the relevant ionic upper levels coupled with 
large photoionization rates (and subsequent cascades). MLD 
postulated that the emission was an artefact of their exclu- 
sion of line blanketing which leads to an overestimate of the 
UV flux and hence of the photoionization rates - further dis- 
cussion of these problems can be found in their paper, while 
their implications for the current dataset are discussed be- 
low. 



sented HST ultraviolet spectroscopy of OB 10-64, suggesting 
approximately B0. 5 lb. From our higher quality optical and 
UV spectra we have adjusted the spectral type sligh tly to 
B0 la (see also Fig. §) . Meanwhile, Massey et al. ( [t986| ) clas- 
sified OB 10-WR1 as a WC6-7 star, which Moffat & Shara 



(1987) revised to WC6+nebL Our high quality spectrum 



3.1 Atmospheric parameters 

A comparison of the spectra of OB 10-64 and HD 167264 in- 
dicates that these stars have very similar, if not identical, 
atmospheric parameters within the errors of our methods. 
This is illustrated in Fig. Fy, where four spectral regions are 
shown for the two stars. These regions contain diagnostics 
of the atmospheric parameters as discussed below: 

Effective temperature: As discussed by MLD, both the 
strength of the Hell features and the relative strengths of 
lines due to Si in and Siiv can be used to estimate the ef- 
fective temperature for spectral types near B0. The former 
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is an excellent diagnostic with a change of effective temper- 
ature of 2000K leading to changes in the Hell equivalent 
widths of a factor of two or three; additionally MLD found 
that the He n line at 4541 A appeared to be well modelled by 
the non-LTE calculation. In Fig. [D, two He II lines at 4541 A 
and 4686A are shown; their strengths seem very similar in 
the two stars. This is confirmed by the observed equivalent 
widths (see below for details of the procedures used to es- 
timates these values) summarized in Table 0. Although the 
line strengths are approximately 20-30% larger for OB 10-64, 
this would correspond to a change in effective temperature 
of less than 500K, such is the sensitivity of the line strength 
to temperature in this T b r range. 

Lines due to Sim (at approximately 4553, 4568 and 
4574A) and Si IV (at 4089 and 4116A) are also shown in 
Fig. hi Again the agreement is excellent as is confirmed by 
the estimates of the equivalent widths listed in Table ti. For 
OB 10-64, the Si ill lines are marginally enhanced and the 
Si IV lines marginally weakened compared with the stan- 
dard star. This would imply that OB 10-64 had a lower effec- 
tive temperature, in contrast to the Hell lines that implied 
that this star was hotter. However, again the differences are 
marginal and would imply a temperature difference between 
the two stars of less than 500K. 

Surface gravity: The profiles of the Balmer lines are 
useful diagnostics of the surface gravity and in Fig. hi, the 
H7 (at 4340A) and H<5 (at 4101A) are shown. In the line 
wings (which have the greatest sensitivity to the gravity), 
the agreement between the two stars is excellent and imply 
that the gravities (assuming that the stars have similar ef- 
fective temperatures) differ by less than 0.2 dex; the major 
uncertainty arises from the difference in S/N of the observa- 
tional data for OB 10-64 compared to the bright standard. 
In the line cores, there are significant differences but this is 
likely to be due to different amounts of rotational or macro- 
turbulent broadening rather than a difference in the sur- 
face atmospheric parameters. Additionally for the OB 10-64 
spectra there were significant amounts of background emis- 
sion from the host galaxy as discussed in the previous sec- 
tion; the discrepancies may in part reflect the difficulty in 
subtracting this background emission accurately. 

Hence we conclude that within the observational uncer- 
tainties, these two stars have the same atmospheric parame- 
ters. MLD deduced an effective temperature of 29000K and 
a logarithmic gravity of 3.1 dex fo r HP 1 67264 from the opti- 
cal spectroscopy of Lennon et al. ( 1992] ; 1993 ) . We have used 
our better quality ISIS observations to check these results 
using the same diagnostics as MLD, namely the Hell lines 
for effective temperature and the Balmer line profiles for the 
gravity. At these relatively high effective temperatures, the 
non-LTE calculations of the Si ill spectrum (and too a lesser 
extent the Siiv spectrum) becomes problematic (see MLD 
for details) and hence these lines were not be considered. 

The observed and theoretical spectra are shown in Fig. 
0, with the latter having been convolved with a Gaussian 
profile with a FWHM of 0.80A to allow for instrumental, ro- 
tational and macroturbulent velocity broadening. In general 
the agreement is good, apart from the Hell line at 4541 A, 
where the wings of the theoretical profile appear to be too 
strong. However even for this case, a decrease of only 1000K 
would lead to a theoretical profile that is significantly weaker 
than the observed. Hence we conclude that the atmospheric 
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Figure 2. Observed (solid line) and non-LTE (dotted line) spec- 
tra for HD 167264. The upper two panels show the He II lines at 
454lA and 4686A which are diagnostics of the effective temper- 
ature. The lower two panels show the gravity sensitive Balmer 
lines, H5 and H7. Excellent agreement is found between the line 
profiles and strengths. 



parameters of MLD are reliable and adopt them for both 
HD 167264 and OB 10-64. Additionally, MLD found that a 
microturbulent velocity of 10 km s _1 was appropriate for 
the Si ill features in B-type supergiant spectra. They also 
noted that the 11 spectra implied a higher microturbulent 
veloci ty - a result found independently by Vrancken et al. 
( 1999J) for early-B-type giants. Here we adopt a microturbu- 
lent velocity of 10 km s _1 for both stars but will consider 
below the effect of a larger value being appropriate. 



3.2 Chemical composition 

The choice of a standard with atmospheric parameters ef- 
fectively identical to those for OB 10-64 simplifies the abun- 
dance analysis. Our procedure was as follows. Firstly metal 
line absorption features were identified in the spectra of 
OB 10-64. For marginal features, the spectrum of HD 167264 
was used as a template and guide. The equivalent widths 
of the features were estimated in both stars by non-linear 
least squares fitting of single or multiple Gaussian absorp- 
tion profiles to the normalised spectra. Normally both the 
positions and widths of the Gaussian profiles were allowed 
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to vary. However for multiplets the relative wavelength sep- 
arations were fixed, while for marginal features in OB 10-64, 
the widths were set to the mean value found for well ob- 
served lines. The equivalent widths are listed in Table Bl 
apart from those for silicon, which are in Table bl all esti- 
mates have been rounded to the nearest 5 mA. 

For most species, these equivalent widths were analysed 
using the non-LTE grids of MLD and the atmospheric pa- 
rameters discussed above. As well as absolute abundance 
estimates a line-by-line differential analysis was also under- 
taken. For the Niii spectrum and one Siiv line, non-LTE 
calculations were not available and LTE methods were used; 
in these cases, the absolute abundance estimates should be 
treated with caution. However, the similarity in both the 
observed line strengths and inferred atmospheric parame- 
ters for OB 10-64 and HD 167264 imply that the differential 
abundance estimates should be reliable. This was confirmed 
by re-analysing the O II line spectrum using an LTE tech- 
nique, where the differential abundances agreed with the 
non-LTE estimates to normally within 0.01 dex. Both the 
absolute and differential abundance estimates are summa- 
rized in Table H, together with the number of features used 
and whether an LTE or non-LTE formulation was adopted. 
Note that the latter information is replicated in Table pi to- 
gether with information on which features were included in 
the abundance analysis. The error estimates in Table ra are 
unweighted standard deviations of the samples; where a sig- 
nificant number of features were measured, the uncertainty 
in the means should be smaller. 
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3.3 Photospheric abundances of OB 10-64 

The atmospheric parameters and chemical composition of 
OB 10-64 appear to be very similar to those of HD 167264. 
Below, we discuss the abundances of individual elements in 
detail: 



Figure 3. Selected spectral regions for OB 10-64 (solid line) and 
HD 167264 (dotted line), showing He I diffuse lines at 4009 and 
4026A, 4471A, 4143A, 4387A. 



Helium: The helium abundance in OB 10-64 has not been 
determined explicitly. This was because most of the He I 
non-diffuse features were either badly blended with metal 
lines (e.g. the line at 4121A) or weak (the lines at 4169A 
and 4437A). The only well observed non-diffuse feature is 
the triplet line at 4713A and as can be seen from Fig. hi 
the profiles for the two stars are very similar. Additionally 
the moderate quality of the observational data for OB 10- 
64 would make reliable profile fitting of the helium diffuse 
lines difficult. However in Fig. H, selected diffuse helium lines 
are shown for the two stellar spectra. In general, agreement 
is excellent and we conclude that within the observational 
uncertainties the helium abundances in the two objects are 
similar. 

Carbon: The carbon abundances are based on a sin- 
gle weak feature (the 4267A doublet) which is of similar 
strength in each star. MLD have commented on the fact that 
the line appears difficult to model properly with the non- 
LTE techniques that we use here. Even in B-type dwarfs, 
calculation of this transition is problematic in both LTE and 
non-LTE. Hence while the absolute abundance estimates 
should be treated with considerable caution, the similarity 
in line strength points to both stars having similar C pho- 
tospheric abundances. 



Nitrogen: This element is particularly important as it pro- 
vides an excellent diagnostic for the mixing of nucleosynthet- 
ically processed material from the stellar core to the surface. 
Hence, a search was made for lines of both N 11 and N ill in 
the spectrum of OB 10-64. Unfortunately all the detections 
were marginal and the corresponding line strengths should 
be treated with caution. The single Nil line implies that 
the two stars have similar nitrogen abundances while the 
four N ill lines implies that OB 10-64 is depleted by approx- 
imately 0.2 dex. However it should be noted that the theo- 
retical N ill line strengths are very sensitive to the adopted 
effective temperature, which coupled with the observational 
uncertainties, implies that this difference is probably not 
significant. 



Oxygen: This is the best observed species with 13 lines 
being analysed. The oxygen abundances in the two stars are 
very similar as can be seen from the absolute abundance 
estimates. It is encouraging that the standard error for the 
differential abundance is significantly smaller implying that 
the use of a differential method is indeed removing some 
sources of systematic error. Assuming that the remaining 
errors are distributed normally, the estimated uncertainty 
in the mean differential abundance would be reduced to ap- 
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Table 5. Absolute and differential abundances for OB 10-64 and HD 167264. Where more than one feature was measured for a given 
species (as denoted by n = number of available lines), the unweighted mean is quoted with the error being the standard deviation of 
the sample. The method of line formation calculation is noted, and for further information see the text. For reference we quote the 
abundances in solar neighbourhood B-type stars from Gies & Lambert (1992) for all elements apart from Mg which is from Kilian (1994) 



ion 


OB 10-64 


§] 

HD 167264 


A[#] 


n 


Method 


Solar B-star 
abundances 


ClI 


7.83 


7.78 


+0.05 


1 


non-LTE 


8.20 


Nil 


7.96 


7.90 


+0.06 


1 


non-LTE 


7.81 


Nin 


8.40±0.07 


8.60±0.27 


-0.20+0.30 


4 


LTE 


7.81 


On 


8.69±0.34 


8.62±0.23 


+0.07+0.15 


13 


non-LTE 


8.68 


Mgn 


7.36 


7.30 


+0.06 


1 


non-LTE 


7.45 


Sim 


7.54±0.18 


7.34±0.12 


+0.20+0.10 


3 


LTE 


7.58 


Si iv 


7.7: 


7.8: 


-0.16 


2 


see text 


7.58 



proximately ±0.02 and hence there is marginal evidence for 
a slight enhancement in OB 10-64. 

Magnesium: The estimates are again based on a single fea- 
ture but the two stars appear to have similar abundances. 

Silicon: The Si ill lines were analysed using both non-LTE 
and LTE methods. The former leads to higher absolute es- 
timates but a similar differential abundance as for the LTE 
calculations. As discussed by MLD, the non-LTE profiles 
of the Si ill lines show some emission in their wings in this 
temperature regime; this is probably an artifact of the lack 
of line blanketing in the calculations. Hence in Table M the 
LTE results are summarized but the absolute abundance 
estimates should be treated with caution. The differential 
abundance implies a small silicon enhancement in OB 10-64 
in contrast to that found from the two Siiv lines. The latter 
result should also be treated with some caution as the anal- 
ysis of one of the Si IV lines was undertaken in LTE (it was 
not included in the non-LTE model ion) , which will make the 
absolute abundance estimates particularly unreliable. Addi- 
tionally the theoretical line strengths are extremely sensi- 
tive to the adopted effective temperature with, for example, 
a change of less than 1000K being sufficient to eliminate 
the discrepancy. Hence we believe that the Si ill differential 
abundance is the more reliable. 

Despite the difficulties for individual ions discussed 
above, we believe that the differential abundance estimates 
should be generally reliable. This reflects the similar at- 
mospheric parameters and metal line equivalent widths for 
the two stars. The former ensures that the method of anal- 
ysis (LTE or non-LTE) is relatively unimportant and in- 
deed we confirmed this by re-analysing the O II equivalent 
widths using LTE methods. Additionally the similarity in 
line strengths makes the choice of the microturbulent veloc- 
ity relatively unimportant provided it is the same in both 
stars; for example, re-analysing the metal line spectra with a 
microturbulent velocity of 15 kms changed the differential 
estimates by less than 0.05 dex. 



3.4 Wind Analysis and Stellar Wind Momentum 

OB 10-64 is an ideal object to test the validity of the con- 
cept of the wind momentum - luminosity relationship. It 




1400 1500 

Wavelength (A) 



Figure 4. A compilation of galactic IUE spectra of two B0 la 
Galactic stars and the STIS spectra of OB 10-64. 



has a well determined chemical composition comparable to 
the solar neighbourhood and an equally well determined ef- 
fective temperature and gravity. The reddening is moder- 
ate and the absolute magnitude follows from the distance 
to M31 (see Table hi). Using a non-LTE model atmosphere 
code with the corresponding effective temperature and grav- 
ity we can calculate the emergent stellar flux and determine 
the stellar radius from the absolute V-band magnitude (see 
Kudritzki 1980). We obtain a radius of R = 34.R©. This 
means that OB 10-64 has stellar parameters very similar to 
the galactic object HD 37128 (e Ori), which has been in- 



vestigated recently by Kudritzki et al. (1999) in their study 
of the wind momentum - luminosity relationship of galactic 
A- and B-supergiants. We therefore expect both objects to 
have comparable wind momenta. HD167264, the Galactic 
comparison object for the abundance study in this paper, 
has a very uncertain distance (and, therefore, radius and lu- 
minosity). From the similarity of the spectra we can only 
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assume that its luminosity must be comparable to OB 10-64 
and HD 37128. 

An inspection of Fig. ]A immediately indicates (albeit 
in a qualitative manner) that the stellar wind properties of 
the three objects cannot be very different. For the quanti- 
tative wind analysis of OB 10-64 we procee d exa ctly in the 
same way as described by Kudritzki et al. (1999). First, we 
measure the terminal velocity v x of the stellar wind from ra- 
diative transfer fits of the strong UV resonance lines of N v , 
C IV and Si IV using the method developed by Haser ( 199E ) 
(see also Haser et al. 1995, Haser et al. 1998, Lamers et al. 
1999). This is similar to the technique applied by Howarth 



et al. ( |1997| ). The fit yields v x = 1650 km s _i for both 
OB 10-64 and HD 167264 (see Fig. |), which is very close to 
the value of 1600 km s" 1 found for HD 37128. 

Subsequently, the mass-loss rate M can be obtained 
from a fit of the observed line profile of Ha. This requires 
a non-LTE unified model atmosphere analysi s incorporat- 
ing the stellar wind and spherical extension (Kudritzki et 
al. 19991) . Assuming a radius of 34 Rq and a terminal ve- 
locity of Voo = 1650 km s" 1 , we calculated a sequence of 
unified models with different mass-loss rates. Fig. M gives 
an idea of the accuracy of the line profile fit, which is some- 
what restricted due to the central contamination by the H II- 
region emission. Despite this, a reasonable determination 
of the mass-loss rate is certainly possible. Values as low as 
l.Ox 10 -6 Moyr -1 can clearly be ruled out together with val- 
ues above 1.75xl0~ 6 M yr _1 . We adopt 1.6xl0~ 6 M©yr _1 
as the model which gives the best compromise fit in Fig. hj 
This means that OB 10-64 has a slightly lower mass-loss rate 
than HD 37128 (2.5 xlO -6 Moyr" 1 , Kudritzki et al. 1999). 
The modified stellar wind momentum log(Mi> 00 (_R/_RQ) ' 5 ), 
which is expected to be a function of luminosity only, is 
28.98 (in cgs-units) compared to 29.15 for HD 37128. Fig. | 
compares the modified wind momenta of OB 10-64 and of all 
t he ga lactic early B-supergiants studied by Kudritzki et al. 
( 1999 ) . Within the intrinsic scatter of the wind momentum 
- luminosity relationship the agreement is excellent. This is 
the first direct demonstration that the concept of the WLR 
works for objects of this spectral type (i.e. early B-types) 
outside our own Galaxy. 



4 WOLF-RAYET ANALYSIS 



4.1 Technique 



We use the non-LTE code of Hillier & Miller (1998) which 



iteratively solves the transfer equation in the co-moving 
frame subject to statistical and radiative equilibria in an 
expanding, spherically symmetric and steady-state atmo- 
sphere. Relative to earlier versions of this code (Hillier 1987, 
1990), two major enhancements have been incorporated, of 
particular relevance to WC-type stars, namely (i) line blan- 
keting, (ii) clumping. 

Our analysis of OB 10-WR1 follows a similar method 
to that recen tly applied to Galactic WC stars by Hillier & 
Miller (1999). Specific details of the (extremely complex) 



model atoms used here for ou r quantitative analysis are pro- 
vided in Dessart et al. (2000), together with information on 
the source of atomic data used for Hel-ll, Cii-iv, Oil- VI, 
Si IV and Feiv-vi. We use an identical (slow) form of the 
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Figure 5. The wind analysis fit (dotted line) for Uoo shown for 
the C IV line only. OB 10-64 is the solid line and HD167264 is 
the dashed. The HD167264 high-resolution IUE spectra has been 
degraded to the STIS resolution. A i>oo=1650 km s — 1 is found. 
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Figure 6. Comparison between the Ha profiles of OB 10—64 (up- 
per) and HD167264 (middle) and HD37128 (lower). The Keck I 
HIRES data of OB 10-64 has been convolved with a Gaussian 
to yield the same resolution as the spectra of HD167264 and 
HD37128 (see Lennon et al. 1993 for details of Galactic star ob- 
servations), and re-binned to a pixel size of 0.3A. The nebular 
lines in OB 10-64 which could not be subtracted properly have 
been snipped out. The profiles of the three stars are qualitatively 
very similar indicating similar wind characteristics. Further, the 
differential photospheric analysis is not compromised by different 
wind properties of OB 10-64 and HD167264. 



velocity law, such that acceleration is modest at small radii, 
but continues to large distances. We assume that the wind 
is clumped with a volume filling factor, /, and that there is 
no inter clump material. Since radiation instabilities are not 
expected to be important in the inner wind we parameterise 
the filling factor so that it approaches unity at small veloci- 
ties. Clumped and non-clumped spectra are very similar, ex- 
cept that line profiles are slightly narrower with weaker elec- 
tron scattering wings in the for mer. Although non-clu mped 
models can be easily rejected (Hillier fc Miller 199S), be- 
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Figure 7. The wind model fit to the Keck I HIRES spectra 
of Ha assuming the parameters outlined in TableW The sharp 
absorption-like features reflect the difficulty or reliably subtract- 
ing the background nebular lines (Ha and [No] 6548.1 A), and 
these parts of the stellar spectra are unrecoverable. Despite this, 
the Ha stellar wind profile can be fitted quite accurately. Three 
NLTE unified model calculations for M = 1.00, 1.50, 1.75 10~ 6 
MQ/yr are shown as dashed curves. 
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Figure 8. Modified Wind Momentum of OB 10-64 (solid) com- 
pared with the galactic supergiants of early B-type spectral type 
in Kudritzki et al. (1999) (open circles). The regression curve is 
a fit to the galactic objects only. 



basis of isolated optical lines from adjacent ionization stages 
of carbon (e.g. Cm A6740/Civ A7700) and/or helium (He I 
A5876/Hell A5412). In practice, this was difficult to achieve 
because of the severe blending in WC winds, so our de- 
rived temperature should be treated as approximate. Cur- 
rent spectroscopic models fail to reproduce the absolute and 
relative strength of the WC classification lines at C ill A5696 
and Civ AA5801-12 (see Hillier & Miller 1998 and Dessart 
et al. 2000). As in other recent spectroscopic studies of WC 
stars, Hell A5412/C iv A5471 were selected as the principal 
C/He diagnostic since the relative strength of these features 
are insensitive to differences of temperature or mass-loss 
rate. Oxygen abundances were impossible to constrai n, since 
the principal diagnostic region spans AA2900-3500 (Hillier 
& Miller 199£). Consequently, we adopt C/0~4 by number 



as predicted by stellar evolutionary models (Schaller et al. 
1992). Based on results obtained for OB 10-64 in Sect, pi we 
shall adopt solar abundances for Si and Fe. 

As discussed by Hillier & Miller (199£), many weak 
spectral features sit upon the continuum so that ideally one 
would wish to compare synthetic models with fluxed spec- 
troscopy. Since absolute fluxes are unavailable for OB 10- 
WR1, we have instead rectified the dataset, taking care to 
use as few 'continuum' regions as possible. 

Our ISIS spectroscopic data of OB 10-WR1 is shown 
in Fig. bl Overall, the spectrum is reasonably well repro- 
duced by our model, except that C in A5696 and C IV 
AA5801-12 are too weak. Nebular emission lines are promi- 
nent in the observed spectrum, and contaminate the Hell 
A6560 and Cm A6740 stellar lines. We find T* ~75 kK, 
logL/L Q =5.68±0.16, v^ ~ 3000 kms" 1 and M/y/J ~ 
lO- 43 M yr _1 , with C/He~0.1 by number, from Hell 
A5412/CIV A5471, despite the low S/N achieved in this re- 
gion from a single 30 minute ISIS exposure. The wind per- 
formance number of OB 10-WR1 is MvooKL/c) ~14 for our 
assumed /=0.1. 

In order to better constrain the stellar parameters 
of OB 10-WR1, one would require UV spectroscopy with 
HST/STIS to permit a reliable reddening determination (the 
chief source of uncertainty in logL/L©), v x via P Cyg reso- 
nance lines, and oxygen abundances from the near-IR diag- 
nostics. For E(B—V)—0.18, the expected UV continuum flux 
from OB 10-WR1 is 2x 10~ 16 erg cm -2 s" 1 A" 1 at 3000A, or 
lxl0~ 15 erg cm" 2 s" 1 A" 1 at 1500A, assuming a st anda rd 
Galactic extinction towards OB 10 (though see Sect. 



cause of the severe line blending in WC winds M /\fj is 
derived by our spectroscopic analysis, rather than M and / 
independently. 



4.2 Spectroscopic analysis of OB 10— WR1 

As usual, a series of models were calculated in which stel- 
lar parameters (T*, logL/Z/Q, Woo M/yff, C/He, O/He) 
were adjusted until the observed ionization balance, line 
strengths, widths, and absolute v— band flux were repro- 
duced. Because of the (substantial) effect that differing 
mass-loss rates, temperatures and elemental abundances 
have on the emergent spectrum, this was an iterative pro- 
cess. 

The wind ionization balance is ideally selected on the 



4.3 Quantitative comparison with HD 92809 
(WC6) 

How do the parameters of OB 10-WR1 compare with 
WCE counterparts in other galaxies? Amongst Galac- 
tic WC6 stars, only HD 92809 (WR23) and HD 76536 
(WR14) a re moderately reddened an d lie at established 
distances ( ]Koesterke fe Hamann 1995 ). Since high quality 
flux calibrated spectroscopy of HD 92809 are available to 
us (Sect. 0), we now carry out an identical analysis of it, 
permitting a quantitative comparison to be carried out. Al- 
though numerous Galac tic W CE stars have been studied by 
Koesterke & Hamann ( |l995| ), including HD 92809, the in- 
clusion of line blanketing and clumping into models has a 
major effect on the stellar luminosity and wind density, jus- 
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Figure 9. Comparison between rectified WHT/ISIS spectroscopy of M31-OB 10— WR1 (solid) and our synthetic model (dotted). Strong 
nebular emission lines severely contaminate the He II A6560 and C III A6740 emission lines, while atmospheric features in the far-red have 
not been removed. 



tifying the need for a re-evaluation of its properties (Dessart 
et al. 2lp00[). 

HD 92809 lies in the Car OBI association, at a distance 
of 2.63 kpc, with a narrow u-band Smith magnitude of 9.71 
mag. Our spectral synthesis of HD 92809, using methods 
identical to those described above, indicates a red dening of 
E(B-V)=0.50 mag (Koesterke & Hamann ( |l995| ) obtained 
E(B-V)=0.38), such that M„ = -4.08 mag, substantially 
lower than OB 10-WR1 (M^ = -5.4 ± 0.4 mag). We adopt 
1-00=2280 km s" 1 f or HD 92809 following the UV analysis by 
Prinja et al. (199C), also somewhat lower than OB 10-WR1. 
The measured emission equivalent widths of selected lines in 
the optical spectra of OB 10-WR1 and HD 92809 are listed 
in Table pi In general, HD 92809 has similar line equivalent 
widths to OB10-WR1, except that the Hen A5412/Civ 
A5471 ratio is much smaller, suggesting a higher carbon 
abundance. Single WC stars appear to show remarkably 
uniform emission equivalent widths at each spectral type, 
oblivious of their stell ar parameters or chem istry. For exam- 
ple, Conti & Massey (Conti & Masseyl989) measured ClV 
A5801-12 equivalent widths in the range 440— 1150A for (ap- 
parently) single Galactic WC6 stars. Therefore, we can have 
confidence that OB 10— WR1 is also single. 

Our spectroscopic analysis was repeated for HD 92809, 
with our final UV and optical synthetic spectrum com- 
pared with observations of HD 92809 in Fig. hd (de-reddened 
by E(B — V)= 0.50 mag following the reddening law of 
Seaton (1979)). Overall, the comparison is excellent for line 
strengths, shapes, as is the match to the continuum distribu- 
tion throughout the spectrum. The number of line features 
that are poorly reproduced is small, and again includes C ill 
A5696, C iv AA5801-12. 

For HD 92089, we derive the following stellar pa- 
rameters: T„=75 kK, logL/L Q =5.3, and M/v7 ~ 



lO" 4 ' 8 Af yr" 1 . Fits to Hell A5412 and Civ A5471 imply 
C/He~0.3 by number, confirming the higher carbon abun- 
dance expected from their line ratio. Near-UV oxygen di- 
agnostics are well matched for HD 92809 with an assumed 
C/Q=4 num ber ratio. As discus sed p reviously by Hillier & 
Miller ([L999J) and Dessart et al. (|2000|) , the inclusion of line 



blanketing into WC models significantly increases stellar lu- 
minosities relative to previous studies - taking the higher 
reddening into consideration, blanketing increases the lumi- 
nosity for HP 9 2809 by ~0.25 dex relative to Koesterke & 
Hamann ( 1995). We confirm the previously derived stellar 
temperature and carbon abundance for HD 92809. A clump- 
ing factor of /=0.1 reduces the global mass-loss rate by ~3, 
so that the overall effect of our study is to decrease the wind 
performance number for HD 92809 to 8 from 108 (Koesterke 
& Hamann 1995). This dramatically alleviates the Wolf- 
Rayet 'momentum problem', such that multiple scattering 
of photons in their highly stratified atmospheres may be 
able to account for their mass-loss properties via radiatively 
driven wind theory. 



5 DISCUSSION 

5.1 Comparison of the Galactic and M31 
B-supergiants 

The differential analysis of the B-types implies that the two 
stars appear to have very similar chemical compositions. 
There is some marginal evidence for a small metal enhance- 
ment in OB 10-64 but this could, for example, just be an 
artifact of an overestimation of its equivalent widths from 
the lower quality spectroscopic data. We believe the differ- 
ential oxygen abundance to be particularly sound given the 
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Figure 10. Comparison between de-reddened (E(B - V)=0.50 mag) IUE/HIRES and MSSSO/DBS spectroscopy of HD 92809 (solid) 
with our synthetic model (dotted). 



13 features used in the analysis and the relatively small 
spread in their results. As these are luminous supergiants 
one must qualify the use of absolute abundance values in 
any external comparison. For oxygen, magnesium and sili- 
con however the abundances derived in HD 167264 (distance 
from Sun ~1.7kpc; Hill, Walker & Yang 1986) are similar to 
the absolute values derived in solar-neighbourhood B-type 
main-sequence stars. The analysis and absolute abundance 
determination of the dwarfs is more reliable and has been 
shown to be in excellent agreement with Hn regions and 
the ISM within ~500 pc of the Sun (Gies & Lambert 1992, 
Rolleston et al. 2000) Hence we can be reasonably confi- 
dent that this agreement means the O, Mg and Si absolute 
abundances of the HD167264 and OB 10-64 can be used. It 
appears that the OB 10 cluster was born from an ISM with 
chemical composition very similar to that of the solar neigh- 
bourhood, and given its proximity to the cen tre o f M31 this 
is a surprising result (see Fig. [13 and Section 5.3). 



The similar nitrogen and carbon abundances in the two 
stars are interesting. MLD in their study of galactic super- 
giants identified three stellar groupings that they designated 
as 'normal', 'processed' and 'highly processed' on the basis 
of their CNO line strengths. By 'processed' they meant stars 
which have C and N line strength ratios which implied that 



their photospheres had been contaminated by the products 
of core CN-cyle burning. HD 167264 was classified as 'pro- 
cessed' as it had strong nitrogen lines in comparison with 
those of carbon. Assuming that the designations of MLD 
are valid, we conclude that OB 10-64 should also be classi- 
fied as 'processed' and hence at least its carbon and nitrogen 
abundances may not be representative of its natal composi- 
tion. 



In Section 2.4 we have assumed a distance to M31 of 
783 kpc and hence have an estimate of Mv = —6.93 ± 0.4. 
The model atmosphere analysis produces a flux distribu- 
tion for the stellar parameters shown in TableW, which pro- 
vides a bolometric correction of —2.9. Hence we determine 
an M bo i = -9.83 ± 0.6, and log L/L = 5.85 ± 0.24. This 
is typical of B-type supergiants (and their p ossib le A-type 
descendants) as found by Kudritzki et al. (1999), and as 
illustrated in FigH. A comparison with evolutionary mod- 
els (see Fig[L3f Schaller et al. 1992) implies an initial main- 
sequence mass for the star of approximately 55M© . Whether 
the star has just come directly from the main-sequence or 
is a blue-loop star with a previous history as a red super- 
giant is unclear. The Schaller et al. models predict that the 
surface CNO abundances should be drastically altered if the 
star has been a red supergiant and is in a blue-loop. While 
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Figure 11. A mosaic of DSS images of M31 tiled together (hence 
the apparent break in the background), with the position of the 
OB 10 association labelled as the square. For comparison, the po- 
sitions of three A-F type supergiants analysed by Venn et al. 
(2000) are shown, with their other star A-207 lying off this field- 
of-view in the outer south west corner of M31. North is up, east 
to the left and the size of the image is 115'x90'. The three crosses 
indicate the positions of the three inner disk H II regions analysed 
by Blair et al. (1982; regions BA75, BA423 and BA289). The 
length of the double arrow shows 5kpc at the distance of M31, 
along the major axis. The projected galactocentric radius of the 
OB 10 association is 5.9 kpc. 



we suggest that it may show some evidence of CNO pro- 
cessing, the quantitative values we derive do not appear to 
support the star having highly processed core material at its 
surface indicative of a red supergiant origin (where the ni- 
trogen abundance is predicted to be approximately +1.0 dex 
higher than normal). Whatever its history, its lifetime is 
likely to have been of the order 3.5±lMyrs (from Schaller 
et al. 1992), in excellent agreement with the age of the more 
evolved WR1 star (with age 3.6 ± 0.9 Myr, see Section 5.2) 
suggesting coeval formation in the OB 10 association. 



5.2 Comparison between Galactic, LMC and M31 
WC stars 

Table M reveals that the stellar and chemical properties of 
OB 10-WR1 and HD 92809 differ substantially, despite orig- 
inating from very similar environments. Following the mass- 
luminosity relatio n for hydrogen-free WR stars of Schaerer 
& Maeder (1992), a present mass of 18T4M0 is implied. 
In contrast, the stellar luminosity of HD 92809 is typical of 
Galactic WCE stars with a present mass of ~11Mq. De- 
pending on the choice of evolutionary model, an age within 
the range 2.7-4.5Myr is implied for OB 10-WR1. 

Since OB 10-WR1 is visually the brightest (apparently 
single) WCE star in M31 (Massey & Johnson 1998), it is 
valid to question whether photometry/ spectroscopy of this 
star is contaminated by bright stars along similar sight lines, 
particularly since 1" corresponds to a physical scale of ~4pc 
at the distance of M31. Excess continuum light would serve 
to dilute the equivalent widths of the WC emission line spec- 



R g (kpc) 



Figure 12. Oxygen abundances as a function of distance from 
the galaxy centre derived in a study of H II regions and the stars 
so far analysed in M31. The open diamonds are the nebular Hn 
region abundances from Blair et al. (1982), and the line shows 
the least-squares-fit to these results, giving a gradient of —0.03 ± 
0.01 dcxkpc~ x . We have recalculated the projected galactocentric 
distances to the Hn regions using an updated distance to M31 
of 783 kpc. The solid circles are the oxygen results from Venn et 
al. (2000), from the three A-F type supergiants shown in Fig. 6. 
The solid star is our photosphcric result from OB 10-64, with the 
error bar representing the standard error in the mean. 



trum, yet from Table pi the emission line strengths of OB 10- 
WR1 agree well with the single WC6 star HD 92809. There- 
fore, in the absence of high spatial resolution imaging we can 
have confidence that the observed continuum flux is princi- 
pally that of the WC star, such that its stellar luminosity 
has not been artificially enhanced. 

Galactic WCE stars are known to span a range of car- 
bon abundances and luminosity - according to Koesterke & 
Hamann (1995) 0.1< C/He <0.5, and 5.0 < logL/L Q < 5.4 
for stars with established distances, and applying a g loba l 
scaling of 0.3 dex to account for blanketing (see Sect. 4.3). 
Therefore, OB 10-WR1 is amongst the least enriched in car- 
bon, with the highest luminosity. These characteristics are 
reminiscent of the nearby WC5+OB binary WR146, re- 
cently analysed by Dessart et al. (2000) who derived log 
L/Lq= 5.7 and C/He=0.08 by number. Extending the com- 
parison to the LMC WCE stars (all WC4 spectral type), re- 
veals quite different conclusions. In contrast with the Galac- 
tic sample, the luminosity of OB 10-WR1 is fairly typi- 



cal of LMC WC4 stars, spanning 5.4 < log L/L e < 6.0, 
with comparable C/He abundances established for Brey 10 
and Brey 50 (see Fig. 4 in Crowther et al. 2000). Therefore, 
OB 10-WR1 is more typical of WCE stars at relatively low 
metallicity. A twice solar metallicity for OB 10 - inferred 
from H II studies (Blair et al. 1982) - would be very difficult 
to reconcile with our derived properties of OB 10-WR1. 
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lOkpc. As noted by Venn et al., neglecting the in- 




log T eff 

Figure 13. The evolutionary tracks from Schaller et al. (1992) are 
plotted for 15— 6OM0 stars, as labeled. The positions of OB 10-64, 
and WR1 are shown. OB 10 appears to have evolved from ~36Mq 
star, whereas WR1 probably evolved from a more massive object. 
Given the existence of an evolved star of lower mass than WR1, 
we can infer limits on the time duration of massive star formation 
in this cluster. 



5.3 Probing the chemical composition and 
evolution of M31 

In Figure ha we have plotted the oxygen Hll region abun- 
dances from Blair et al. (1982) showing the abundance gra- 
dient they derive, and the position of OB 10-64. We have 
recalculated the de-projected galactocentric distances of all 
the objects in this plot by assuming a distance of 783 kpc, an 
inclination angle of 77.5° , major axis angle of 37.5° , and cen- 
tre of M31 to be RA(2000) = 00 h 42 m 44"4, Dec. (2000) = 
41° 16' 08" 97 (the latter from Hjellming & Smarr 1982). The 
galactocentric distance of OB 10 is 5.9 kpc, and the Blair et 
al. results suggest that the gas metallicity at this position 
should be significantly above solar (~9.1dex). Our abso- 
lute abundance of 8.7 dex is well below this value, and while 
we acknowledge the difficulty in comparing absolute abun- 
dances from different methods we are confident that OB 10- 
64 does not have a metallicity significantly higher than this 
value — the star has metal-line strengths and atmospheric 
parameters very similar to those of t he Galactic standard 
HD 167264 as discussed in Section ! 



5.1 



Venn et al. (200C) have presented an analysis of three 



other stars in M31 (A-F supergiants) and the oxygen abun- 
dances for these are also plotted in Figure ft% Fig. [Ill also 
shows the spatial positions of the A-type supergiants and 
OB10 in the M31 disk. We see that the present results 
strengthen the suggestion of Venn et al. that the abun- 
dance gradient of oxygen is quite flat, extending the stel- 
lar results into the supposedly metal rich region interior to 



nermost Hll regions results in good agreement between the 
nebular and stellar results, albeit with a negligible gradient, 
so it is important to try to understand the current discrep- 
ancy in the inner region. In the H n region work the observed 
quantity used to derive [O/H] is R23, the ratio QOII]A3727 + 
[OIII]AA4959, 5007)/H/3, which is observed to increase with 
R 9 . This ratio has been calibrated as a function of [O/H] 
by a number of authors, and indeed as Venn et al . pointed 
out, using the calibr ation of Zaritsky et al. ( 1994 ) instead 
of Pagel et al. (1979), reduces the [O/H] abundance gradi- 



ent from Blair et al.'s value of —0.03 dex kpc - to —0.02 
dex kpc -1 . Indeed, within the range of R23 defined by the 
Hll regions in question, 0.0 < R23 < 0.8, there are signif- 
icant differences between the slopes and absolute values of 
the various R23 versus [O/H] relationships (see Figure 11 
of McGaugh 1991). In fact McGaugh's own calibration re- 
sult s in a n even flatter gradient than that of Zaritsky et 
al. (1994), and systematically lower oxygen abundances by 
0.1 to 0.2 dex, with the largest difference being at higher 
metallicities (see Kobulnicky et al. 1999, equations 9 and 
10 for the parameterization). For metallicities above solar 
there are almost no H 11 regions for which there are indepen- 
dent estimates of the electron temperature due to the weak- 
ness of the [OIII] A4363 line. Thus the various calibrations 
in use depend either upon extrapolated empirical relation- 
ships, or upon models for which one must make certain as- 
sumptions (for example concerning abundances, abundance 
ratios, model atmosphere fluxes and dust). 



Fortunately, Galarza et al. (1999) have presented new 
results for M31 H 11 regions between 5-15 kpc, and three 
of these are very close to the OB association OB 10 (re- 
gions designated K310, K314 and K315). K315 is spatially 
coincident with OB 10 and we assume that it is sampling 
the nebu la excited by its OB stars. Interestingly, Galarza et 
al. (1999) discuss their Hll regions by morphological type; 
center-bright, ring, complex and diffuse, but only those re- 
gions defined as center-bright exhibit a significant gradient 
in i?23- Of the regions mentioned above, only K314 is a 
center-bright region, it is about 1.5 arc minutes away from 
OB 10-64, and McGaugh's calibration gives 12+log[0/H] = 
8.98, whereas Zaritsky 's gives 9.17. We can apply these cal- 
ibrations to all three Hll regions mentioned above, plus the 
5 sub- apertures centered on K315 for which Galarza et al. 
also give data, which results in oxygen abundances 



( h_999| ) 



ranging from approximately 8.7 to 9.0 dex using McGaugh's 
calibration, or from 8.9 to 9.2 dex using Zaritzsky's (see 
Table H). The former calibration appears to give absolute 
results in better agreement with the stellar results, indeed 
our absolute abundance of 8.69±0.1 agrees within the errors 
with the McGaugh calibration of 8.7-9.0 dex. We therefore 
use this relationship to rederive oxygen abundances using 
the published Hll data from Blair et al. (1982), Dennefeld 
& Kunth (1981) and Galarza et al. (center-bright regions 
only). A discussion of which calibration is more physically 
reliable is outside the scope of this paper, however it seems 
consistent to homogeneously re-derive the abundances from 
the full observational dataset using a single relation. We have 
chosen the McGaugh calibration simply because it provides 
better agreement with the OB 10 stellar results, but would 
caution against using any one of these calibrations blindly 
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Figure 14. All Hn region results from Blair et al. (1982; dia- 
monds), Dennefeld & Kunth (1981; triangles) and Galarza et al. 
(1999; crosses) using the calibration of McGaugh (1991). Stel- 
lar points are represented by filled symbols as in Figure 8. The 
straight line represents a least squares fit to the nebular data 
giving a gradient of —0.018 ± 0.01 dexkpc. 



before more work is carried out on their reliability and the 
reasons for discrepancies. 

These results are shown in Figure u4 together with the 
stellar results and the least squares fit to the nebular results 
which gives an oxygen abundance gradient of— 0.018L0.01 
dexkpc -1 , with an intercept at R g = 0.0 of 9.02 dex. Note 
that the corresponding values using Zaritsky's calibration 
are —0.025 dexkpc -1 and 9.22 dex, however the gradient de- 
termined from Hn regions in the inner, possibly metal rich 
regions, still depends critically upon the functional form of 
the calibration used. Clearly this problem can only be re- 
solved by more detailed observations and modelling of these 
inner regions. The stellar results imply that the [O/H] gra- 
dient is extremely flat in the range 5 to 25 kpc. We have not 
fitted a gradient to the four stars as they are not a homoge- 
neous dataset in terms of their atmospheric parameters and 
analysis methods; the star analysed here is a B-type super- 
giant, those at 12 kpc are A-type supergiants, while that at 
22 kpc is an F-type supergiant. Certainly, more stellar ob- 
servations are clearly required, as any gradient which does 
exist may be masked by statistical scatter or real variations 
at any particular R 9 . 

Without the analysis presented here for OB 10-64 we 
would have assumed that the initial metallicity of OB 10- 
WR1 was much higher than solar, and that the wind analysis 
and stellar evolution discussion would have been based on 
quite different parameters. This indicates that one must be 
careful in assuming an initial metallicity determined from 
global galactic properties when characterising the evolution 
of single evolved stars, especially WR stars. It is necessary 
to use more direct methods to determine the state of the 
natal gas in the environment of the cluster formation. 



6 CONCLUSIONS 

We have carried out detailed, quantitative analyses of high- 
quality spectral data of two massive stars in the young OB 10 
association in M31. This represents the most detailed study 
of hot, massive stars in this galaxy to date, and shows that 
quantitative extragalactic stellar astrophysics is not only 
possible with these objects, but is a powerful diagnostic tool 
advancing our knowledge in several related areas. 

Our analyses allows us to probe massive stellar evolu- 
tion, mass loss and stellar wind characteristics in the most 
massive, luminous stars in the Local Universe. We show that 
we can determine accurate chemical abundances in B-type 
supergiants which are probes of the current metallicity of 
the ISM. Our stellar analysis indicates that the empirical 
calibration of the R23 ratio in Hn regions at metallicities 
of solar and above is not well constrained. Hence when dis- 
cussing stellar evolution from characteristics of evolved mas- 
sive stars (WR stars, LBV's, blue and red supergiants) one 
should be careful to use direct metallicity determinations of 
the starforming region - photospheric abundance estimates 
in massive B-type (and A-type, see Venn et al. 2000) super- 
giants provide an excellent means of doing so. 

We have demonstrated that it is possible to analyse the 
physical and chemical properties of individual WC stars in 
M31 using 4m class telescopes, although more precise lumi- 
nosities and abundances await UV spectroscopy. Compar- 
ison between OB 10-WR1 and WCE stars in other galax- 
ies, particularly HD 92809 in our Galaxy, reveals that its 
high luminosity and low C/He ratio is more typical of WCE 
stars in metal poor environments, in spite of the low num- 
ber statistics. Since OB 10 is in the inner region of the M31 
disk (R 9 = 5.9 kpc), we would have wrongly concluded that 
WR1 originated in a very metal rich region, had we assumed 
a metallicity from the previous Hn region analysis. 

We have carried out a detailed wind-analysis of the B- 
type supergiant and (assuming previously determined dis- 
tances to M31) find its wind momentum-luminosity rela- 
tionship to be in excellent agreement with similar Galactic 
early B-type stars. This is a further indication (cf Kudritzki 
et al. 1999), that a properly calibrated wind momentum- 
luminosity relationship within the Local Group may allow 
accurate distance moduli to be determined to the most lumi- 
nous stars in galaxies within ~20 Mpc. In further papers we 
will analyse a larger set of M31 A and B-type supergiants, 
extending the radial baseline from the central to regions to 
R s k, 20 kpc. 
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Table 3. Equivalent widths, W\, for metal lines in OB 10-64 and 
HD 167264. Typical errors for OB 10-64 are 10% for lines with 
W-a > 70mA and 20% for lines with W\ < 70mA. Measurements 
of W\ for HD 167264 should normally be accurate to 5% 

Wavelength Species W\ (mA) Status 

OB 10-64 HD 167264 



Table 3. continued 
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On 


4277.40 


On 


4277.90 


On 



400 

90 

115 
115 

80 
120 



4317.14 


On 


75 


4319.63 


On 


85 


4319.93 


On 




4349.43 


On 


115 


4366.89 


On 


135 


4379.05 


Nm 


125 


4414.90 


On 


110 


4416.97 


On 


135 


4481.13 


Mgn 


105 


4481.33 


Mgn 




4510.92 


Nm 


55 


4514.89 


Nm 


85 


4523.60 


Nm 


10 


4590.97 


On 


115 


4595.96 


On 


100 


4596.17 


On 





35(1 



75 



110 



Blended 



non-LTE 


non-LTE 


Blended 



non-LTE 



non-LTE 



70 


non-LTE 


80 


non-LTE 


150 


non-LTE 


LOO 


non-LTE 


L10 


LTE 


135 


non-LTE 


60 


non-LTE 


95 


non-LTE 


90 


LTE 


L20 


LTE 


45 


LTE 


90 


non-LTE 


70 


non-LTE 



Wavelength Species Star Number Status 

OB 10-64 HD 167264 



4609.44 


On 


4610.14 


On 


4630.54 


Nil 


4631.38 


Si iv 


4640.64 


Nm 


4641.82 


On 


4641.90 


Nm 


4647.42 


Cm 


4649.14 


On 


4650.25 


Cm 


4650.84 


On 


4651.47 


Cm 



4661.63 

4673.74 
4676.24 



On 

On 
On 



115 
60 

230 

520 

875 

75 
170 



70 
55 
295 

105 
810 



non-LTE 
non-LTE 
Blended 

Blended 

Blended 



85 non-LTE 

135 non-LTE 



Table 4. Stellar atmospheric and wind parameters derived for 
OB 10-64. 



Parameter 


Value 






T e g 


29000 


± 1000 


(K) 


logg 


3.1 


± 0.2 


(dex) 


« 


10 


± 5 


(km s -1 ) 


M 


1.6xl0~ 6 


±°;§ X 10-6 


(Meyr 1 ) 


Voa 


1650 


±70 


(km s -1 ) 


R* 


34 


±6 


(R©) 


logL/Lg 


5.85 


±0.24 


(dex) 



Table 6. Emission equivalent widths (W\ in A) for selected lines 
in the optical spectra of OB 10-WR1 and HD 92809, with typi- 
cal errors of ±10%. More uncertain values for OB 10— WR1 are 
shown in parenthesis when strong nebular emission severely con- 
taminates the stellar spectrum 



Spectral 


feature 


W x (A) 




Major 


Minor 


OB 10-WR1 HD 92809 


C III AA4647-50 


He II A4686 


813 


881 


Hen A5412 




53 


30 


C iv A5471 




32 


10 


Ov AA5572-5604 


Om A5590 


35 


45 


C m A5696 




145 


198 


Civ AA5801-12 


He I A5876 


1085 


907 


Hen A6560 


Civ A6560 


(130) 


150 


C III AA6727-60 




(100) 


215 


C iv A7063 


O iv AA7032-53 


135 


140 


C iv A7724 


C iv A7737 


165 


190 
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Table 7. Comparison between stellar parameters of OB 10-WR1 
in M31 with the Galactic WC6 star HD 92809 (WR23). Assuming 
M„ is appropriate, typical uncertainties are ±10kK, ±0.12 dex, 
±100 km s — 1 , ±0.05 dex, ±15% respectively for each parameter. 



Star T, 
kK 


logL 


kms 


log M 

MQyr" 1 


C/He 


mag 


OB10-WR1 75 
HD 92809 75 


5.7 
5.3 


3000 
2280 


-4.3 
-4.8 


0.10 
0.30 


-5.4 
-4.1 



Table 8. Hn regions observed by Galarza et al. (1999) in the 
vicinity of OB 10; note that K315 is coincident with the position 
of OB 10-64 for which we obtain an oxygen abundance of 8.7 dex. 
Derived abundances [O/H]^ and [O/HJjvf refer to the use of the 
calibrations of Zaritsky et al. (1994) and McGaugh (1991) respec- 
tively. 



H II region 


Morpholgical 


log(i?23) 


[0/H]z 


[0/H] M 


identifier 


type 








K310 


diffuse 


0.47 


9.00 


8.85 


K314 


center-bright 


0.24 


9.17 


8 


98 


K315 


complex 


0.39 


9.07 


8 


90 


K315a 


diffuse 


< 0.23 


> 9.16 


> 8 


98 


K315b 


diffuse 


0.23 


9.16 


8 


98 


K315c 


diffuse 


0.16 


9.20 


9 


01 


K315d 


ring 


0.49 


9.00 


8 


83 


K315e 


ring 


0.63 


8.84 


8 


69 
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